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\D . Abstract 
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Anisotropic flow of kaons and antikaons is studied in heavy-ion collisions at CERN 
SPS and BNL RHIC energies within the microscopic quark-gluon string model. In 
the midrapidity range the directed flow of kaons v\ differs considerably from that 
of antikaons at SPS energy (E^t, = 160 AGeV), while at RHIC energy (y/s = 130 
AGeV) the excitation functions of both, kaon and antikaon, flows coincide within 
the statistical error bars. The change is attributed to formation of dense meson- 
dominated matter at RHIC, where the differences in interaction cross-sections of 



o 



o 

kaons and antikaons become unimportant. The time evolution of the kaon anisotropic 
flow is also investigated. The elliptic flow of these hadrons is found to develop at 
midrapidity at times 3 < t < 10 fm/c, which is much larger than the nuclear pass- 
ing time t pass = 0.12 fm/c. As a function of transverse momentum the elliptic flow 
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increases almost linearly with rising pt- It stops to rise at pt > 1.5 GeV/c reaching 



the saturation value v^ipt) ~ 10% 



0. 
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1 Introduction 



The main aim of experiments on heavy-ion collisions at relativistic and ultra- 
relativistic energies is to study the properties of nuclear matter under extreme 
conditions and to search for distinct signals from the deconfined quark-gluon 
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plasma (QGP) (see [1] and references therein). The transverse collective flow 
of particles is at present one of the most intensively studied characteristics of 
heavy-ion collisions [2,3,4], because the flow is directly linked to the equation 
of state (EOS) of the system. If even a small amount of the QGP is formed in 
the course of the collision, it would lead to a reduction of pressure [5,6,7] and 
a softening of the EOS [8,9] that can be detected experimentally. To study the 
properties of transverse particle flow the method of Fourier series expansion 
[10,11] has been proved to be very useful: 
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Here pt, y, and <j) are the transverse momentum, rapidity, and the azimuthal 
angle of a particle, respectively. The unity in square brackets represents the 
isotropic radial flow, while the other terms are refer to anisotropic flow. The 
first Fourier coefficient in Eq. (1) v\ = (cos0) is called directed flow. It repre- 
sents the averaged ratio of the particle momentum along the impact parameter 
axis to the transverse momentum, v\ = (p x /pt)- The second Fourier coefficient 
i>2 = (cos (20)) is called elliptic flow. It characterizes the eccentricity of the 
ellipsoid of the particle azimuthal distribution, v 2 = ((p x /Pt) 2 — iPy/Vt) 2 )- 
Methods proposed for measurement of azimuthal anisotropies in heavy ion 
collisions can be found in Refs. [12,13,14]. 

The anisotropic flow is a function of rapidity, transverse momentum, and the 
impact parameter of an event b, i.e. v n = v n (xj), where {0^=1,2,3} = {y,Pt,b}. 
Therefore, the following differential distributions are usually applied 
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The idea that the elliptic flow can carry important information about the early 
stage of heavy-ion collisions has been discussed already in Ref. [15]. This sug- 
gestion is supported by macroscopic hydrodynamic and microscopic transport 
simulations, which show that elliptic flow saturates quite early [16,17,18,19,20], 
while directed flow develops almost until the stage of final interactions [21,22]. 
However, the directed flow of hadrons with high transverse momentum can be 
used as a probe of hot and dense phase of the collision [23] due to the early 
freeze-out times of these particles. 

In the present Letter we are studying the anisotropic flow of kaons, K +y s and 
K 0, s, and antikaons, K~'s and K 's, produced in lead-lead (CERN SPS) and 
gold-gold (BNL RHIC) collisions at E lab = 160 AGeV and v^ = 130 AGeV, 
respectively. Our investigation is inspired by several reasons: First of all, in 



heavy-ion collisions at GSI SIS (1 - 2 GeV) and BNL AGS (2-11 AGeV) 
energies the excitation functions of the directed K + and K~ flow are found to 
be different. At SIS energies the K + flow is zero or slightly anti-correlated with 
the nucleon flow [24,25]. This effect becomes even more pronounced at AGS 
where a strong anti-correlation between the K + and the nucleon flow appears 
[26]. In both cases, this anti-correlation can be explained by the influence 
of an in- medium kaon potentials derived from effective chiral models [27]. 
Kaons experience a weak repulsive potential as a result of superposition of a 
repulsive vector potential and an attractive scalar potential, whereas antikaons 
experience a strong attractive potential. Therefore, the directed K + flow at 
these energies is anti-correlated (antiflow) [28,29] while the K~ flow is expected 
to be similar to the nucleon flow (normal flow) [30]. 

Does it mean that at energies of SPS and higher, where the mean fields do 
not play a dominant role, the anisotropic kaon flow would be similar to that 
of antikaons? Also, multiple particle production at high energies is usually 
described in the Monte Carlo microscopic models via the formation and break- 
up of strings stretching between quarks and diquarks (or antiquarks). Kaons, 
that can be formed on the so-called leading quarks (u and d quarks belonging 
to primary hadrons), are K + (us) and K° (ds) only, but not K~ (us) and 
K (ds). This may affect the flow of i^'s, because leading hadrons are carrying 
larger transverse momenta. Can we see the difference between the coefficients 
vi£ +K and v% +K in different rapidity intervals? Finally, K~'s and K 's can 
be absorbed via the channels such as K~ + p — > A + 7r + , etc., whereas there 
are no analogous reactions for K +, s and K 0, s. How important is this reaction 
asymmetry at SPS and, especially, at RHIC energies, where the matter is 
expected to be meson-dominated? 

To answer these questions we employed the microscopic cascade quark-gluon 
string model (QGSM) [31,32,33,34]. Description of the model is given in Sec. 2. 
Rapidity and transverse momentum dependences of directed (Sec. 3) and ellip- 
tic (Sec. 4) flow of kaons and antikaons are calculated in Pb+Pb and Au+Au 
collisions at SPS and RHIC energies, respectively. Time development of the 
both flows is studied as well. Finally, conclusions are drawn in Sec. 5. 



2 The model 



Similar to the dual parton model (DPM) [35], VENUS [36], and the string 
fusion model (SFM) [37], QGSM is based on Gribov-Regge theory (GRT) [38] 
accomplished by a string phenomenology of particle production in inelastic 
hadron-hadron (hh) collisions. To describe hadron- nucleus and nucleus- nucleus 
collisions the cascade procedure of multiple secondary interactions of hadrons 
was implemented. The model incorporates the string fragmentation, formation 



of resonances, and rescattering of hadrons, but simplifies the nuclear effects 
neglecting, e.g., the mean fields or evaporation from spectators. As indepen- 
dent degrees of freedom QGSM includes octet and decuplet baryons, octet and 
nonet vector and pseudoscalar mesons, and their antiparticles. The momenta 
and positions of nucleons inside the nuclei are generated in accordance with 
the Fermi momentum distribution and the Woods-Saxon density distribution, 
respectively. Two particles can interact inelastically with the cross section 
&inei( s ) or elastically with the cross section a e i(s) with the "black disk" prob- 
ability P(s, b) = 9 (a(s)/n — b 2 ). Here a(s) = <Ji ne i(s) J i-a e i(s) is the total cross 
section, b is the impact parameter, and 9 is the step function. Pauli blocking 
of occupied final states is taken into account. 

Several different subprocesses are responsible for the string formation in in- 
elastic collisions. A single string can be produced as a result of the valence 
(di) quark annihilation. This subprocess is important in hh interactions at 
low energies since the annihilation cross section drops as s -1 / 2 . The baryon- 
antibaryon annihilation in QGSM is modeled via the string junction annihi- 
lation with the creation of three strings. In other subprocesses two or more 
strings can be produced as a result of the colour exchange mechanism, while in 
single or double diffractive scatterings one or two strings are formed due to mo- 
mentum transfer. The Pomeron, which is a pole with an intercept «p(0) > 1 
in the GRT, corresponds to the cylinder-type diagrams. The s-channel dis- 
continuities of the diagrams, representing the exchange by n-Pomerons, are 
related to process of 2k (k < n) string production. If the contributions of all 
n-Pomeron exchanges to the forward elastic scattering amplitude are known, 
the AGK cutting rules [39] enable one to determine the cross sections for 2k- 
strings. The hard gluon-gluon scattering and semi-hard processes with quark 
and gluon interactions are also incorporated in the model [34]. The inclusive 
spectra in the QGSM automatically have the correct triple-Regge limit for 
Feynman variable x — > 1, double- Regge limit for x — *• 0, and satisfy to all 
conservation laws [31]. 

Because of the uncertainty principle hadrons produced in string fragmentation 
can interact further only after a certain formation time. However, hadrons 
containing the valence quarks of the primary nucleons can interact promptly 
with reduced interaction cross section. To calculate the statistical weights for 
the different subprocesses the experimental total, elastic, and inelastic cross 
sections are used. If the experimental data are lacking, the additive quark 
model relations, the one-pion exchange model, detailed balance considerations, 
and isotopic invariance are employed. 

The variety of subprocesses included in the model, - diagrams with quark 
annihilation and rearrangement, associated with the exchange of Reggeons, 
subprocesses with colour exchange which are connected to (multi) Pomeron 
exchanges, as well as hard processes with large momentum transfer, - is quite 



rich. The contribution of each subprocess to the inelastic cross section depends 
strongly on the bombarding energy. This circumstance enabled the model to 
predict quantitatively the directed flow of nucleons at SPS [5] more than six 
years before the first measurements [40,41] and to describe the characteristics 
of the elliptic flow of charged particles at RHIC [42] measured quite recently 
[43,44,45]. The formation of a QGP is not assumed in the present version of the 
QGSM, although one may consider the strongly interacting field of coloured 
strings as a precursor of the plasma. Our strategy is to find clear quantitative 
and qualitative deviations of the experimental results from model predictions. 
These deviations can then be taken as an indication of new physics, most 
probably the creation of QGP, not included in the model. 



3 Directed flow 



Time evolution of directed flow of kaons and antikaons in minimum bias 
Pb+Pb collisions at 160AGeV is presented in Fig. 1. Here the coefficient vf(y) 
is calculated in different transverse momentum intervals at early, t = 3 fm/c 
and t — 10 fm/c, and the final stage of the reaction, t > 60 fm/c. To avoid 
ambiguities, all resonances in the scenario with early freeze-out were allowed 
to decay according to their branching ratios. The flow evolution is seen quite 
distinctly. At early stages of the collision directed flow of both kaons and 
antikaons is oriented in the direction of normal flow similar to that of nucle- 
ons [23]. Within the error bars there is no differences between (K + + K°)'s 
and (K~ + K )'s. At this stage the matter is quite dense, mean free paths 
of particles are short, and similarities in kaon production and rescattering 
dominate over inequalities caused by different interaction cross-sections. It is 
worth mentioning that the directed flow of K 's and ICs is already sizeable 
at t — 3 fm/c. This can be explained by a kick-off effect associated with the 
early stage of the collision, when the nuclei pass through each other. Later 
on the system becomes more dilute. For both kaons and antikaons the di- 
rected flow experiences significant transformations as seen in Fig. 1. Already 
at t = 10 fm/c the antiflow of antikaons starts to built up in the midrapidity 
range. Note, that t — 10 fm/c corresponds to the maximum of the kaon dN/dt 
distribution over their last elastic or inelastic interaction [46]. Here the differ- 
ences in interaction cross sections and possible reaction mechanisms become 
crucial. As was found in previous investigations [47,48], particles emitted in 
the direction of dense baryon-rich matter (i.e., normal flow) will interact more 
frequently compared to those emitted in the opposite, i.e. antiflow, direction. 
This circumstance reduces the resulting directed flow of kaons almost to zero 
in a broad midrapidity range. Due to larger interaction cross-sections of i^~'s 



and K 's with other hadrons, the directed flow of these particles changes the 
orientation from a weak normal to strong antiflow. Even (K~ + K )'s with 



high transverse momentum demonstrate distinct antiflow, while the flow of 
(K + + K°ys remains almost unchanged compared to that at t — 10 fm/c. 

The directed flow of kaons and antikaons in minimum bias Au+Au collisions 
at yt~s = 130 AGeV is shown in Fig. 2 again at early stages, t = 3 fm/c and 
t = 10 fm/c, and at the final one, t > 100 fm/c. It is interesting that at t — 3 
fm/c (i) the flow of (K + + K°)'s coincides within the statistical errors with the 
(K~ + K ) flow, and (ii) the flow is generally very similar to that at the SPS 
energy at time t = 3 fm/c. Except of the target and projectile fragmentation 
region, where again the flow is probably produced by the initial kick, the kaon 
flow at this early stage of gold-gold collisions at RHIC energy is isotropic 
with respect to the impact parameter axis. Note that the colliding nuclei have 
passed through each other already at times less than 0.12 fm/c. The spatial 
anisotropy in the distribution of baryonic charge seems to be unimportant at 
this stage. At t = 10 fm/c not only the directed flow of antikaons, but also 
that of kaons becomes antiflow-aligned at midrapidity. Similar behaviour has 
been found within the RQMD model for the directed flow of nucleons at RHIC 
[49], suggesting that the nucleon directed flow is a side effect of the elliptic 
flow. The flow of produced particles, pions [49,50] and kaons [50], was found 
to be very flat at \y\ < 2, in stark contrast to the QGSM predictions. We are 
awaiting the experimental data to resolve this problem. 

At the final stage of the reaction the excitation functions of both kaon and 
antikaon flow have a similar antiflow behaviour in the midrapidity range com- 
pared to the zero flow of (K + + K°)'s at lower energies. However, the kaon flow 
at high transverse momenta 0.9 < pt < 2.0 GeV/c is elongated in a normal di- 
rection. Its slope at the stage of thermal freeze-out is steeper than the slope at 
t = 10 fm/c, i.e., the flow of high p r kaons becomes stronger. This means that 
at RHIC energies the nuclear medium is rather dense, the net baryon density 
is small, and that kaon-meson (including resonances) reactions dominate over 
kaon-baryon interactions, which can significantly modify the flow of K~'s and 
K 's at SPS energies. 



4 Elliptic flow 



The study of the development of elliptic flow of kaons in minimum bias Pb+Pb 
collisions at 160 AGeV, displayed in Fig. 3, reveals that there is no particular 
direction of kaon emission in the transverse plane at t — 3 fm/c. From this 
time on until the freeze-out stage one cannot distinguish between the elliptic 
flow of kaons and antikaons in the various p r intervals. The resulting ellip- 
tic flow of .ff's and K's is positive, i.e., kaons are propagating preferentially 
in the reaction plane. The flow increases with rising transverse momentum. 
This behaviour appears to be quite similar to those obtained from the hy- 



drodynamic calculations. Recall that in the hydro models the strong in-plane 
elliptic flow of hadrons is explained exclusively by the so called almond shape 
of the overlapping region of the colliding nuclei. As a consequence of this initial 
spatial anisotropy, in-plane pressure gradients should be stronger than their 
out-of-plane counterparts [15]. The fluid acceleration is proportional to the 
pressure gradient divided by total energy density. During the early stage of 
the expansion this ratio remains roughly constant. When the spherical shape 
of the expanding region is restored, the pressure becomes uniform. As was 
argued in [15], the typical time scale for the transverse flow is the trans- 
verse size of the system divided by the sound velocity in the medium. For a 
semi-peripheral gold-gold or lead-lead collisions this criterion gives us approx- 
imately 7 fm/c. At early times the elliptic flow in hydrodynamics increases 
linearly with time. Therefore, it is quite plausible that the elliptic flow of 
(anti)kaons at midrapidity in the QGSM calculations is formed within the 
time interval 3 fm/c < t < 10 fm/c in the center-of-mass (CM) system. 

Similar features are observed in the development of the elliptic flow of kaons 
and antikaons in Au+Au collisions at y/s = 130 AGeV, as one can see in Fig. 4. 
Their elliptic flow is close to zero in the midrapidity range \y\ < 2 at t = 3 
fm/c. Even the flow of high-p t particles, which are produced in very energetic 
hadronic collisions at the beginning of the reaction and which decouple from 
the system earlier than other particles, shows the same "isotropic" tendency. 
However, the elliptic flow of kaons and antikaons in the midrapidity range 
2 < y < 4 (SPS) and \y\ < 1.5 (RHIC) seems to be formed between 3 fm/c 
and 10 fm/c, i.e., the flow indeed probes the early stage of heavy-ion collisions. 
The final elliptic flow of both, (K + + K°)'s and (K~ + K )'s, is positive in 
accordance with the predictions of Ref. [15] and its strength increases with 
rising p t . 

The dip at y — 0, clearly seen in Fig. 4, looks a bit peculiar. However, the 
signal is not unique. The appearance of similar midrapidity dip in the i>2(y) 
distribution of hadrons, calculated in [50] at full RHIC energy ^/s = 200 AGeV 
within the UrQMD model, is linked to the particle formation time. Our study 
supports this conclusion. Figure 4 shows that the effect is stronger for high-pr 
kaons, which have longer formation time due to higher momentum. Shorter 
formation times lead to increase of the elliptic flow, and the dip vanishes. 

The transverse momentum dependence of the elliptic flow in both reactions 
is presented in Fig. 5. Again, in the very early phase of the collision the 

flow v 2 yPt) is absent. At time t = 10 fm/c the elliptic flow of kaons and 
antikaons seems to be already formed. It is close to zero for kaons and antikaons 
with low transverse momenta, p t < 0.25 GeV/c, and then rises linearly up to 
sfipt) ~ 10% within the interval 0.25 < p t < 1.5 GeV/c. At transverse 
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momenta higher than 1.5 GeV/c the elliptic flow saturates in accord with 
experimental results [51]. Note, that if the hydrodynamic regime would be 



reached, the excitation function of the elliptic flow v^ipt) should grow further 
almost linearly above p t > 1.5 GeV/c also [17]. 

The delay in the development of the elliptic flow and the saturation of the flow 
after a certain p r limit can be explained as follows. The bulk production of 
new particles takes place within the first few fm/c's, when parton and hadron 
collisions are extremely energetic. Kaons, as well as other hadrons produced 
as a result of string fragmentation, are emitted isotropically in the azimuthal 
plane. To develop the anisotropic flow these hadrons have to rescatter in a 
spatially anisotropic nuclear matter. But because of the uncertainty principle 
the only particles allowed to interact immediately (although with reduced 
cross sections) are those containing the valence quarks. For all other hadrons 
a proper formation time should pass between the rescatterings. Due to the 
finiteness of colliding nuclei, there is a saturation of number of secondary 
interactions per hadron with the transverse momentum higher than certain p t . 
Therefore, the elliptic flow of high-p t particles also saturates. Here it is worth 
to mention that the problem of the formation time of hadrons produced in a 
very dense and hot nuclear medium, which is formed in heavy ion collisions 
at energies of RH1C and higher, is studied insufficiently A shorter formation 
time will reduce the mean free path of particles, thus increasing the collision 
rates and enhancing the elliptic flow. This can be the reason why the plateau 
in the excitation function v 2 (Pt), predicted by the QGSM (see Fig. 5), is 
a bit lower compared to the data. This question should be clarified in nearest 
future. 



5 Discussion and conclusions 



In summary, the anisotropic flow of kaons and antikaons is studied in lead- 
lead and gold-gold collisions at SPS and RHIC energies, respectively, within 
the microscopic cascade quark-gluon string model. The directed flow of kaons 
v 1 + (y, Ap t ) at SPS energy is close to zero in a broad rapidity range except 
of the projectile and target fragmentation regions, while the directed flow of 
antikaons vf +K (y, Ap t ) has a linear antiflow slope. At early times t < 3 
fm/c for both particle species the directed flow is oriented in the direction of 
normal flow. The difference in final distributions is attributed to the different 
interaction cross sections with the surrounding nuclear matter and to the large 
absorption cross section of antikaons with baryons. In heavy-ion collisions at 
RHIC energies a dense meson-dominated matter with low net baryon density is 
formed. Hence, the directed flow of kaons becomes similar to that of antikaons, 
i.e., in both cases the flow is elongated in antiflow direction at \y\ < 2 and in 
normal flow direction at higher rapidity, \y\ > 2. The predicted antiflow of i^'s 
in the central rapidity range is an order of magnitude stronger compared to 



the signal predicted by the UrQMD model. These are the main results of our 

study. Since the differences between the functions vf +K (y) and v^ +K (y) 
in the fragmentation regions in different transverse momentum intervals are 
found to be within the statistical error bars, we conclude that the formation of 
fast kaons on leading u and d quarks plays only a minor role at RHIC energies. 

The elliptic flow of both kaons and antikaons in the midrapidity range is built 
up at times between 3 fm/c and 10 fm/c, although the passing time even 
for central heavy-ion collisions is 1.44 fm/c and 0.12 fm/c at SPS and RHIC 
energy, respectively. In contrast to the midrapidity range, the elliptic flow of 
kaons and antikaons with rapidity y > 1.5 continues to develop after t = 10 
fm/c in the center-of-mass system. The rapidity distribution of the function 

v 2 is flat over a broad range \y\ < 3 at RHIC energies. As a function 
of transverse momentum the elliptic flow of K's and K's rises linearly with 
increasing p t and saturates at p t > 1.5 GeV/c. The plateau at high transverse 
momenta can be explained by saturation of number of secondary interactions 
experienced by hadrons with a transverse momentum higher than a certain 
limiting value. 
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Fig. 1. The directed flow v% of kaons (solid circles) and antikaons (open circles) in 
minimum bias Pb+Pb collisions at SPS energy as a function of rapidity y in high 
Pt interval 0.6 < pt < 0.9 GeV/c (left panels) and for all transverse momenta (right 
panels) at times t = 3 fm/c (upper row), 10 fm/c (middle row), and final (bottom 
row). 
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Fig. 2. The same as Fig. 1 but for minimum bias Au+Au collisions at RHIC 
(s 1 ' 2 = 130 AGeV). For high j^-kaons transverse momentum varies from 0.9 GeV/c 
to 2.0 GeV/c. 
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Fig. 3. The same as Fig. 1 but for the elliptic flow vi of kaons (solid circles) and 
antikaons (open circles) in minimum bias Pb+Pb collisions at SPS. 
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Fig. 4. The same as Fig. 2 but for the elliptic flow vi of kaons (solid circles) and 
antikaons (open circles) in minimum bias Au+Au collisions at RHIC. 
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Fig. 5. Transverse momentum dependence of the elliptic flow of kaons (solid circles) 
and antikaons (open circles) in minimum bias Pb+Pb collisions at SPS (left panels) 
and Au+Au collisions at RHIC (right panels) at early times t = 3 fm/c (upper row), 
10 fm/c (middle row), and freeze-out times (bottom row) of the reaction. 
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